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Abstract 


Nowadays, the behaviour of fluids in weightlessness is a great interest for a lot of scien- 
tists all around the world. A particular aim of this study is to predict how fluids behave 
in space. The use of cryogenic propellant - liquid hydrogen ( LH 2 ) and liquid oxygen 
(. LO 2 ) - is increasing, especially in space launcher. The main questions the designers of 
these launchers have to solve is how heat transfer are taking place in absence of gravity 
and how does the gravity affect the boiling crisis. 

Experiments have been performed at the Low Temperature Laboratory (SBT) of the 
CEA/GRENOBLE in 2009 to answer this question for liquid oxygen. The experiments 
that are carried out in this master thesis are based on the same principle. However, 
the studied fluid is LH 2 . It consists on measuring the thermal transfers in nucleate and 
him boiling regime for different levels of gravity, until Og. The compensation of gravity 
is done by using a huge magnetic held, leading to an important gradient of this held, 
used to levitate the hydrogen. Such a facility is available at the CNRS of Grenoble in 
the National Laboratory of the Intense Magnetic Fields (LNCMI). This resistive coil 
of 20MW delivers a magnetic held of 15 Teslas. 

For various thermal heat hux, the temperature difference between the wall and the 
huid has been measured to plot the boiling curve : Nukiyama curve. These experi- 
ments have been carried out for several level of gravity and thermodynamics conditions 
(pressure and temperatures have been studied independently). 

Results were very satisfying and the main goal of this study has been achieved: these 
results are able to answer the main questions of designers for the new Vinci engine, 
used in the next generation of spacecraft Ariane. 
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1 Introduction 


1.1 Context 

Studies on fluid are of a great interest for many scientific teams. They enable to predict 
the behaviour of a fluid and then to size as precisely as possible each component of 
a system. The following study is focused on cryogenic fluids used in space launchers: 
liquid hydrogen (L// 2 ) and oxygen (L0 2 ) [1], 

Thermal behaviour of such fluids are already well known on earth, which enabled, 
to build space launchers used in different space programs. However, some designers 
wish to know how these fluids behave in absence of gravity, to be able to use the 
launchers in space after being put off for a while, when the pipes are hot. 

To answer this question, it is of pair importance to perform experiments at zero gravity. 
Several ways of getting the compensation of gravity can be used, as described in the 
background part below. 

1.2 The study 

In this master thesis, Nukiyama curves [2] are established to define boiling behaviours 
for different thermodynamics condition of LH 2 and different level of gravity. These 
curves compare the evolution of the wall temperature due to an imposed flux. This 
study was made by using the magnetic compensation solution regarding the gravity 
problem. 

To fulfil the needs, a specific cell has been built with an internal heater whose purpose 
is to transfer the desired power directly to the fluid. An observation system is installed 
to see, record and understand what happens. In the cell, a cooling circuit regulates the 
temperature around 20K, corresponding to the liquid hydrogen temperature at desired 
pressure (1 and 2 bars). It turns out that obtained results are very promising. This 
report present all the work I’ve done during this master thesis. 

The first section of this report is composed of a background part, for the reader to 
understand the main aspects of this study. The next section is composed a presen- 
tation of different materials and methods that have been used. Results of this study 
are then presented before coming to the discussion part where several parameters are 
presented, and finally the last part concludes this master thesis. 


1 



2 Background 

In this chapter, several aspects of the problem are presented, in order for the reader 
to be familiar with the main concepts of the thesis. First of all different techniques 
of getting microgravity are explained. A description of the magnetic compensation 
method is then presented, and finally an overview of thermal heat transfer in fluids is 
described 

2.1 Classical techniques to obtain microgravity 

Several techniques are nowadays utilized. All of them having their pros and cons, here 
is a brief review of these technique. 

2.1.1 Drop tower 

The drop tower uses an inertial technique to compensate earth gravity. A freely falling 
body experiences gravity forces and an opposing inertial force, appearing due to the 
change of reference frame: from a Galilean frame attached to earth’s surface to a non- 
Galilean frame attached to the body accelerating at the acceleration due to gravity. 
The result is that the body experiences zero g. 

In drop tower, vacuum conditions are established, in order to reduce the aerodynamic 
forces. Experimental installations are put in shuttles falling down in the tower. By 
using this technique, it is possible to obtain a gravity of the order of 10 -5 g during 4 
to 5 seconds. 



Figure 1: Drop tower 


It’s difficult to increase this time, since on the one hand the height of the tower would 
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increase as the square of microgravity time and on the other hand, the resulting kinetic 
energy would be so large that a very huge brake system is required to dissipate this 
amount of energy . Moreover, the rotation of the earth itself should also be taken into 
account. 

2.1.2 Parabolic flights 

CNES and ESA have a modified plane called Airbus Og, which enables to carry out some 
experiments in microgravity. This method enables to compensate gravity to the order 
of 10 _2 g during 20 seconds. A typical parabolic flight consists of about thirty parabolic 
trajectories. However, between two phases the plane has to go back to the previous el- 
evation, during which an over gravity of 1.8g is applied on the plane. Moreover, such a 
technique is rather expensive and cannot be applied to LO 2 and LH 2 for safety reasons. 



Figure 2: Airbus Og 


2.1.3 Space shuttle 

Experiments could be conducted on space shuttle board when in orbit. In an orbit 
around earth on a reference frame fixed to the space shuttle, the centrifugal force 
cancels out gravity forces thus establishing zero g conditions. Typically gravity levels 
of the order of 10 -4 can be achieved. 

2.1.4 Space station 

Space stations, such as MIR erstwhile and the International Space Station enable to 
obtain a microgravity of the order of 10 -2 to 10 -3 g. The centrifugal force due to the 
rotation of the station around earth in an orbit is exactly in opposition to its weight. 
Residual acceleration due to human activities, vibration and constant repositioning to 
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Figure 3: Space shuttle 


be kept in orbit induce oscillations (called jitter) and thus affect the quality of zero g 
obtained. However, the duration is much longer : a lab can be kept in microgravity 
during the whole life of the station (10 to 20 years). 

2.1.5 Magnetic compensation 

Zero g can be achieved by using strong magnetic fields through a technique called 
magnetic compensation. In this method the magnetic properties of the material to 
be levitated (diamagnetism /paramagnetism) are exploited to exert a strong magnetic 
body force in opposition to the gravity force. When sufficiently strong, the gravity 
forces can be cancelled to achieve zero g. 

This is a very useful method to obtain microgravity since: 

• The installation is on earth and immobile, exerting no space / mass constraints on 
the system, and could be easily secured against hazard risks related to L0 2 and 
LH 2 

• Microgravity can be theoretically obtained for an infinite time duration. 

• Several level of gravity can be achieved, to simulate for example the gravity 
conditions of other planets (such as Mars, Moon, etc...) 

• Time variations of gravity are possible 

• This technique is available anytime. 

• The experiments can be less expensive than the other techniques of zero g, de- 
pending solely on the current in the coil (which can be considerably reduced by 
using superconducting coils) 

This technique has however some drawbacks : volumes that could be levitated are 
small and the strong magnetic field may have some influence on the fluid one want to 
study. 


4 


2.2 Gravity compensated by a magnetic field 

In this section, the basic principles of magnetic compensation are explained. 


2.2.1 Principle 

In order to compensate earth gravity for diamagnetic fluids, a magnetic held is applied. 
The interaction between this held and the huid leads to a force per unit volume given 
by the following expression: 

~f = B 2 (1) 

2/i 0 

where y represents the magnetic susceptibility of the huid , /j 0 the permeability of the 
vacuum and if the magnetic held, where is the gradient operator. According to the 
magnetism convention, y is dimensionless. However, the mass susceptibility y p or the 
molar one xm is sometimes used in the literature. M is the molar mass. 

The relation between these 3 values is : 

X = PX P = j^Xm (2) 

The studied huid is Hydrogen, which is a diamagnetic material. This kind of material 
is characterised by a very small, and negative susceptibility y. For liquid Hydrogen, 
it corresponds to (— 10 -6 ) whereas for its vapour phase it goes from (— 10 -9 ) to (— 10 -6 ). 

Eq (1) shows that this force exists only if 

• there is a magnetic held 

• this held varies in space 

A volume element considered in a magnetic held is subject to two forces: 

• its weight : = p~$ 

• the magnetic force : ~E\ n = p^-X^ B 2 

2p, 0 

By applying the fundamental principle of dynamics, it leads to 

F^ + t g = pt*pf !L XB 2 = ^ (3) 

to simplify expressions, let’s call the compensation vector B 2 . Hence, by dehning 

the vertical unit vector such as = — \g\v% and .vt z = G z the exact magnetic 
compensation is obtained for : 

(4) 

Xp 


Three remarks can be made: 



— y This expression shows that for diamagnetic material (x < 0), Go is oriented in 
the same direction as the gravity vector. 

-4 To levitate a diamagnetic material such as hydrogen, huge field’s gradient have 
to be applied. For comparison, a paramagnetic material such as oxygen, has a 
susceptibility 10 3 times higher than hydrogen, thus requires a smaller gradient to 
compensate gravity for the same volume 

-4 the exact compensation for liquid Hydrogen at 20K is obtained for 
G 0 = -991.47’ 2 .m-' 


2.2.2 Inhomogeneity of the compensation field 

It is useful to define a new gravitational vector representing the effective compensation: 




( 5 ) 


If the magnetic force is not applied gt = if, which corresponds to the normal earth 
conditions. If the magnetic force is equal and opposite to the gravitational force, the 
residual gravity is null: the microgravity conditions are fulfilled. 

The residual gravity can be considered as the effective gravity the fluid is subjected to. 
This method is adapted to compensate exactly tiny volumes. Indeed, the exact com- 
pensation is realised only at two specific points [3]. Since the gravitational force is 
constant at the scale of the installation, but the magnetic held is not, it leads to some 
non homogeneities of compensation: each point being subject to a slightly different 
gravitational vector. 


A simple solenoid can create only two points where the gravity is exactly compen- 
sated. However, only one of these point is stable -i.e: if a small volume spreads from 
this point, it will automatically come back to its initial position, whereas for an unstable 
point it will tend to move further away (see explanation page 9). 

To characterise the lack of compensation in the volume, it’s useful to define the non 
homogeneity vector at each point: 


Cr — Go 


This vector represents the residual forces applied to the fluid. 
Mixing Eq (4),Eq (5)and Eq (6) leads to 


( 6 ) 


gt 

~W\ 


( 7 ) 


which can be analysed as the non dimensional residual gravity vector. 
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2.2.3 Magnetic compensation in a simple solenoid 

2.2.4 Vector analysis 

With a simple solenoid, the geometry is axisymmetric and the induction vector if and 
if are defined only with r, in the radial direction of the coil and z components, in the 
axis direction. Hence, for our problem, one gets in cylindrical coordinates: 


= ^B 2 = 




It leads to the following expression for e : 

— > 1 . 1 G r (r,z) 

G 0 | G z (r, z) — G 0 

Moreover, on the z axis, it leads to some interesting properties due to symmetry: 




Go V dz G ° 


Since a non- zero gradient is necessary to get the spatial compensation, zero g is possible 
only at the two extreme positions of the coil, where |Cr z | is maximum as it can be seen 
in figure 4: These calculations have been carried out with the Levitation software 
developed by the High Magnetic Field National Laboratory (LNCMI) of the CNRS 
Grenoble. 

For a specific material, a theoretical expression describes the relation between the 
magnetic held, the inhomogeneity level and the volume that can be levitated: 


D 1 / 3G 0 R 


(8) 


This expression is given by considering the second order Taylor expansion of B in 
the expression of and using Maxwell equations : div ( if )=0 and rot( if) = if and 
considering a vertical magnetic held.R corresponds to the characteristic length of the 
levitated volume (in m) 


For example, to levitate hydrogen, a magnetic held of 5.5T can compensate gravity 
with a residual level of 1% for a radius of 1.22 mm, which corresponds to a volume of 
7.6 mm 3 of hydrogen 

The line drawn in hgure 4 represents the compensation G z for hydrogen. It can be 
easily seen that only 2 points can be compensated. Since the magnetic held depends 
directly on the current howing through the coil, three conhguration can be presented: 
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z (m) 

Figure 4: Magnetic induction B z and Compensation G z on z axe of the coil for Hydro- 
gen 


• if the current is not high enough, the \G Z \ curve doesn’t cross the compensation 
line: there is no compensation at all 

• if the current increases until the \G Z \ curve becomes tangent to the compensation 
line, there is exactly one compensation point 

• if the current increases even more, there are two compensation points 

The third configuration is the only one which is interesting to compensate the gravity 
in a volume. The part where the curve is above the compensation line is called the 
overcompensation zone, whereas when the curve below this curve is called the under 
compensation zone. It can be noticed that the green line corresponds to the opposite 
of the gradient. That’s why it’s inverted compared to figure 5 

2.2.5 Axial stability 

The major difference between a paramagnetic and a diamagnetic fluid is the sign of the 
vector G due to the contribution of the magnetic susceptibility. Hence, to compensate 
the gravity in a paramagnetic fluid, Cr has to be positive in order to get a magnetic 
force upwardly oriented (if the gravity is down oriented) whereas for diamagnetic fluid, 
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such as hydrogen, G has to be negative since Xm is negative too. 

There are only two forces acting on a bubble: upthrust buoyancy and the magnetic 
force. For each fluid near its compensated area, two cases are possible as it can be seen 
in figure 5: 



• The bubble is at C: a small upward displacement (right side on the z-axis of 
figure 5) puts the bubble in the overcompensated area, the gravity becomes lower 
leading to a smaller upthrust buoyancy. The resulting force tends to move the 
bubble down, back to C. If the small displacement is directed in the other way, 
the bubble would be in the under-compensated zone, where the gravity becomes 
higher, leading to a stronger upthrust buoyancy. The resulting force tends to 
move the bubble back to C. This point is stable for vapour since whatever the 
displacement, the bubble come back to it’s initial position. 

• The bubble is at D: a small upward displacement put the bubble in the under 
compensated area, where the upthrust buoyancy is stronger than the magnetic 
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force, the resulting force tends to move the bubble even upper. If the small 
displacement is directed in the other way, the bubble would be in the overcom- 
pensation zone, where the gravity is lower, leading to a smaller buoyancy. The 
resulting force tends to move the bubble even lower. This point is unstable for 
the vapour phase. 

It can be noticed that the unstable point for vapour phase D is stable for liquid 
phase and the stable point for vapour C phase is unstable for liquid phase: for 
vapour the magnetic force is directed downward to compensate the buoyancy (the 
weight of vapour is negligible) and for liquid it’s directed upward to compensate 
the weight 

2.2.6 Radial instability 

according to (7) 


That is why the study of G r enables to determine the radial force applied on a bubble. 
A simple taylor expansion at first order shows that e r is proportional to r. Hence, if 
a small displacement deviates the bubble from the r-axis, we get: G r > 0, leading to 
e r < 0 since G 0 < 0 and finally g* < 0 

Hence radial acceleration g* is oriented to decreasing r, thus the resulting upthrust 
buoyancy is oriented to increasing r. That’s why hydrogen is radially unstable 


9% = 9 | ^ 


thus 


G r 


(9) 



L 



exact comp 
^ movement 

movement 


exact compensation point 


movement of the liquid 


movement of the vapor 


Figure 6: radial instability of an hydrogen bubble 


2.2.7 Modification of the gravity level 

Lots of combinations of current (I) and axial position (z) are possible. Let (io, zq) be 
one such combinations. 
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By projecting Eq(5 ) on z axis, it leads to: 




( 10 ) 


dB 2 

where G z (I,z) — — — (/, z) and g* is the residual gravity, g the earth gravity, xp the 
mass susceptibility of hydrogen, and g 0 the magnetic permeability of vacuum and G z 
the axial component of the magnetic compensation force. 

By definition, we get r/*(/ 0 , zo) = 0 , hence 


9 — G z (Iq, zq) 
Zno 


The presence of a metallic insert could have change the shape of the magnetic held. 
However, in our case no metallic insert is added and thus the magnetic held is directly 
proportional to the current in the solenoid, that is why 

with f and function independent of the current. Hence it leads to 


9 = 9*{I,z) +x P I 2 J(z) 

and 


( 11 ) 


9 = X P Io 2 -f(zo) 

If now the altitude z 0 is kept constant , the relation becomes 

g = g*(I,z 0 ) + x P I 2 -f(z 0 ) 

Since 

4 - =*,/(*»), 

1 o 


g = g*(I,z 0 ) + I 2 - 


9o 


It finally gives the relation between the current and the residual gravity : 

g*(rzo) _ j _ 

9 h 


( 12 ) 


T = 'l i 

I 0 




(13) 


This relation will be useful to know the applied current corresponding to the desired 
level of gravity. 
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2.3 Heat transfer in boiling fluids 

2.3.1 Three modes of heat transfer 

There are three fundamental mechanisms of heat transfer : conduction, convection and 
radiation. 

Conduction 


This mode corresponds to a transfer of heat in a material without any physical move- 
ment of material. Fourier’s law gives the heat flux density q CO nd i n W.nT 2 transmitted 
by this phenomenon: 

(Icorul — —X^T (14) 

where A (in W.rrf 1 . K' 1 ) is the thermal conductivity of the material and T (in K) is the 
temperature held. 

Convection 


For this mode, the heat transfer is due to a movement of particles. In the general 
case, the heat hux density transmitted by convection from a solid bulk to a liquid is 
given by 

q CO nl = h(T b - Tf)lt (15) 

where h (in W.rrf 2 .K~ l ) is the convection exchange coefficient, T p (in K) is the temper- 
ature of the solid wall, Tf is the temperature of the fluid far away from the solid and 
it the normal vector to the solid directed outward. 

The transfer coefficient h depends strongly on the configuration of the system. Two 
different types of convection are distinguished: the natural convection, where the sys- 
tem is put in a force held, such as gravitational force, or forced convection where the 
huid is put in motion through external forcing via another system, a pump for example. 

Radiation 


The heat hux density q ra d emitted by a system at a certain temperature T is given 
by 


<jfd = e.a B -T , -Tt 


(16) 


with e the emissivity of the system, a B = 5.6703 * 10~ 8 W.rrf 2 .K~ 4 the Stefan-Boltmann 
constant, T (in K) the temperature of the system and it the normal vector of the 
system directed radially outward 
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2.3.2 Boiling 


Boiling happens when the ambient pressure is larger than the vapour pressure. It is a 
very effective way of heat transfer. Two boiling type exist [4]: 

• The homogeneous one, where the fluid is boiling in all the liquid. This type is 
very difficult to achieve experimentally 

• The heterogeneous one, where the fluid is boiling in contact to a hot wall, in a 
different material than the studied fluid. 

The experiment that have been carried out in this master thesis have been done with 
this second type. Extreme care has been taken to ensure that : 

• The fluid is pure 

• The liquid is in equilibrium with its vapour or undercooled. 

• The volume is closed 

• The heating heat flux was given by a solid material 


There are two different forms of boiling that will be further explained in the next 
part : nucleate/pool boiling, and film boiling. 

Pool boiling 


To study pool boiling, two different types of experiment are carried out, presented 
in figure 7: 

• using a hot wire [5] 

• using a heating flat plate 

Our experiment have been carried out with the flat plate mode. The heating part 
is a hot wall placed in the bottom of the cell in which the fluid is filled. In normal 
condition, due to gravity, the liquid is at the bottom part of the cell while the vapour 
is above it. The heating surface is then totally in contact with the liquid. 

Dividing the power injected in the liquid by the area of the heating plate gives the 
heating flux density. 
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Figure 7: experimental devices : hot wire (on the left) and plat plate (on the right) 
T sat is the saturation temperature of the fluid 


2.3.3 Nukiyama curve 

Several parameters can be measured: the heat flux density q (in W.m ‘ 2 ), the saturated 
temperature of the fluid T sat and the temperature of the wall, T p . The heat flux density 
q is imposed and the temperatures are measured. Each experimental point is recorded 
in stationary regime, when all parameters are stables. The overheating is defined as 
AT sat =T P — T sat . The curve q=f(AT sat ), called Nukiyama curve can then be plotted 
(see figure 8) 



Figure 8: Nukiyama curve 


14 


This curve starts when the fluid is at saturated temperature (A T sat = 0), where 
the imposed heat flux is null. 

From A to B 


In the first region, there is no vaporization. The temperature of the fluid in contact 
with the heater is below the boiling point of the fluid. There is only conduction and 
natural convection in the fluid. The heat flux is transmitted directly to the liquid very 
slowly by a conduction process first and then by a convective one. from B to C it’s the 
onset of nucleate boiling. 

From C to D 


As soon as the overheating achieves the threshold at B, that is when the tempera- 
ture of the first layer of liquid in contact with the heater is equal to the boiling point, 
bubbles appear. As it can be shown in the graph, the temperature of the wall strongly 
decreases when the nucleate boiling starts (B-C), due to the latent heat of vaporiza- 
tion: Since this phase change is endothermic, it absorbs heat while creating bubbles. 
This boiling continues until reaching the point D. Even if the heat flux constantly 
increases, the overheating is quite small (i.e, the wall stays quite cold) , due to the 
efficiency of this thermal transfer. As the heat flux is increased more and more bubbles 
are created (figure 9). 



From D to E 


As soon as the threshold D is achieved, the heat flux is too large to be evacuated 
and for the nucleate boiling to keep going . The corresponding flux is called CHF 
(critical heat flux). There is so much vapour that the liquid can’t reach the wall to 
cool it: a vapour film constantly covers the wall (figure 10). 

Thermal transfer is highly impacted by this vapour layer as it has low thermal conduc- 
tivity. The temperature of the wall (A T sat ) quickly increases an reaches a very different 
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equilibrium value, even if the heat flux stays constant. 



Figure 10: Film boiling. 


From E to F 


It is film boiling. When the heat flux is still increased, the wall temperature keeps 
increasing with a quite high temperature. The slope of the curve is very small, mean- 
ing that the thermal transfer is very poor: a small increase in the flux causes a huge 
increase in the wall temperature. If the heat flux is sufficiently high, the material can 
reach its melting point and then destroy the wall (in F). 

There is an hysteresis: 


• If the heat flux decreases during the nucleate boiling phase, experimental points 
stay in the nucleate boiling part of the curve 

• If the heat flux decreases during the him boiling phase, experimental points don’t 
come back to the nucleate boiling curve but stay in the him boiling one (F-G) 

There is a minimal heat hux that can be reached in him boiling (G). At this point, 
the liquid starts to wet the wall again and then the nucleate boiling phase can be re- 
established. 

From F to G 


This part of the curve is very difficult to obtain due to the instability of the change 
between the two modes, and cannot be achieved with this method. To get it, the 
overheating has to be imposed , not the heat hux. 

This boiling curve is universal. Whatever the studied compound, the shape of this 
curve is the same. This thesis will be focused on the plotting of this curve for liquid 
hydrogen for different pressures and levels of gravity. 
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2.3.4 Nucleate boiling 

To start boiling, an energy threshold has to be crossed. The temperature of the heater 
brings the needed amount of energy to do it. However, several other parameters also 
play important role in it, especially surface defects, which can considerably reduce the 
threshold. That’s why the needed energy to start boiling on a polished surface is higher 
than on a rough one [4], 

Each defect is a potential nucleation site, where bubbles will be able to develop. It’s 
a geometrical imperfection (hole), containing vapour or air. While the threshold is 
crossed, this site become active and bubbles start to nucleate from it. Different shape 
of site leads to different behaviour of boiling [6] . 

There are several steps in the creating of bubbles: 

• waiting 

• inertial growing 

• thermal growing 

• leaving 

the following figure show this different steps. 


mmm 

Waiting 



Inertial growing 



Thermal growing 





Leaving 


Figure 11: developement of a bubble 


Waiting 


A small volume of vapour stays in the nucleation site. The liquid phase, at saturated 
temperature, recovers this site, replacing the area occupied by the previous bubble. 
The heating compound heats up the liquid by conduction. During this step, the bub- 
ble doesn’t grow. 

Inertial growing 


As soon as the thermal limit layer is built back, the overheating becomes sufficient 
for a new bubble to start growing. Pressure forces are then higher than the surface 
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tension, leading to the fast development of the bubble with an hemispherical shape, 
the micro layer between the bubble and the wall enables to improve the heat transfer 
by evaporating the liquid caught in it. 

the bubble grows according to the Rayleigh-Plesset model [7] with a characteristic 
speed: 



(17) 


where r b is the radius of the bubble, A P = P v — P L is the pressure difference between 
the liquid and gas phase, p L is the liquid density, and t the time. 


Thermal growing 


As soon as pressure forces are of the same order of magnitude as that of the surface 
tension, the bubble gets a spherical shape, to minimize its energy. Its growing is due 
only to the vaporization of the liquid at the interface. Heat transfer is then done only 
by diffusion at the interface and grows as [7], even if a small flow is present: 



where Al is the thermal conductivity of the liquid phase, AT is the overheating, p v 
is the liquid density, H is the latent heat of vaporization, and Dt,l is the thermal 
diffusivity of the liquid phase. This expression is correct only in a hot liquid, far away 
from the wall. 

Leaving 


As soon as the bubble is as long as the capillary length, 


r b ~ 



(19) 


where a is the surface tension, g the gravity and A p the difference of density between the 
liquid and gas, the bubble leaves due to upthrust buoyancy. It brings some saturated 
liquid with it , leading to a perturbation in the thermal boundary layer. A new cycle 
can then start again. 


2.3.5 Critical Heat Flux 

As explained above, the heat transfer is more intense in nucleate boiling than in him 
boiling. Moreover, the hotter the wall, the better the transfer. Hence, it’s very im- 
portant for an industrial heat transfer device to work on C-D part of Nukiyama curve, 
close to point D. 
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However, this point shouldn’t be crossed, otherwise it would heat up the wall due to 
low conductivity of the the vapour layer. The transition occuring at CHF is, by defi- 
nition, critical: it’s instantaneous, unstable, chaotic and irreversible (even if the heat 
flux is decreased directly after having crossed this point). 

A lot of experiments have been done since the discovery of the CHF, trying to un- 
derstand this phenomenon. However, even if lots of configurations have been studied, 
industrials are dependent only on empirical data. Hence they have to apply a quite 
important security coefficient to size their installations. Knowing exactly what hap- 
pens would reduce this coefficient and then costs. 

This background section has presented a brief overview of concepts that are neces- 
sary to understand and treat the subject of this master thesis 
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2.4 Information about hydrogen 

Hydrogen gas can be considered as an ideal gas. It is a colourless, odourless, tasteless, 
non-toxic, non-corrosive, non-metallic diatomic gas at standard pressure and temper- 
ature. Its density is very low, which makes it necessary for any practical applications 
to either compress the hydrogen or liquefy it. Figure (12) shows its phase diagram 



Figure 12: Phase diagram of hydrogen 


As shown in the diagram, at 1 and 2 bar, the corresponding phase change is around 
20K. That is why cryogenic installation is necessary to carry out our experiments. 

It is an highly diffusive and highly buoyant gas, which mixes with air very quickly 
upon release. Hydrogen coexists in two different forms, ortho and para hydrogen, the 
fraction of each of these components depends on the temperature. Normal hydrogen 
at room temperature is 75% ortho (nuclear spins aligned) and 25% para (spins anti- 
aligned). At 20 K, there are 99.821 % of para and 0.179 %ortho . It takes 3-4 days 
until a new equilibrium state is reached. 

However, magnetic impurities and also small oxygen concentrations are able to catalyse 
ortho-para conversions raising the rate by several orders of magnitude to the order of 
hours. Catalysts such as Fe(OH ) 3 or even small concentration of oxygen can accelerate 
this transformation and/or enable to get a desired fraction of each. 

There are only few differences between physical properties of the two spin states. The 
most important one is their large energy difference which results in major differences 
for the specific heats and thermal conductivities. 
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3 Methods and materials 


In this section, the whole installation is described. First, an explanation about different 
functions of each component is given, and then the sizing of some of them is presented. 

3.1 The experimental setup 

The installation is called LHYLA (Liquid HYdrogen Levitation Apparatus). A sketch 
of it is presented below on figure 13. It shows the different parts of the installation, 
which are precisely described later on. 



to the 
outside of 
the building 


Figure 13: Global drawing of the installation 


It can be noticed that the vision system is not represented on this figure, but is 
presented later. 
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3.1.1 The cryostat 


The creation of undesired ice is the main enemy of cryogenists. Indeed it could break 
down a whole installation. Whatever the substance, if its solidification temperature is 
higher than the colder part of the installation, it could create ice. 

That is the reason why cryogenics is always associated with vacuum technique: it pre- 
vents from having undesired elements, and also removes convection movements in the 
cryostat which would impact the cooling by adding this thermal transfer to others. 

To cool down an installation, the heat has to be extracted. It can be done either by 
direct or indirect technique. The first one consists in directly extracting the heat of 
the device by thermodynamic cycles, and the second one consists in using cryogenic 
fluids such as liquid helium or nitrogen for example, to extract heat from the device. 

The first method is much more expensive than the second one. That’s why for our 
purpose, the second one has been chosen. 

The cryostat is the compound which enables to thermally isolate the low temperature 
device (fluid and mechanical parts) to the external atmosphere. 

Our cryostat is represented in figure 14: 



Figure 14: Moving of the cryostat with an overhead crane 
As shown it’s a quite large structure (3m high) which can be manipulated only 
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with overhead crane. Such a set up is designed to study both ergol components of any 
rocket engine, i.e. L0 2 and LH 2 . As explained before, the microgravity zone for LH 2 
is situated at the top site of the levitation coil; while in the case of L0 2 , this zone is 
at the bottom of the coil. 



Figure 15: Inside of the cryostat (on the left) and top of the cryostat(on the right) 

To limit heat exchanges inside the cryostat, the three modes already presented are 
utilized: 

• Convection: A secondary vacuum (1CT 5 mbar) is maintained inside it by a turbo- 
molecular pump. This prevents from any gas exchange. 

• Conduction: Thanks to the secondary vacuum, heat losses by conduction are 
only possible through metallic pipes going from top of the cell (endoscope, gas 
feeding). However, they are very long and thin and then the heat transfer is 
insignificant (it will be precisely expressed in the sizing part). 
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• Radiation: a thermal shield (not represented on the previous figure) totally sur- 
rounds the experimental cell. It is cooled by gaseous helium and achieves a 
temperature around 40 K. Hence it reduces a lot radiating heat transfer (from 
5W to some mW) 

To see what happens inside the cell, a vision system is installed, as shown in figure 
(15). This system is further described. The flange situated on top of the cryostat (see 
figure 15) contains: 


• 2 electrical motors for the rotation and translation of light and vision mirrors 

• 2 electrical motors for the translation of light and vision endoscope 

• Electrical connection of sensors and heating resistance 

• Cooling pipes: cold helium at 6K enters into the cryogenic path system to regulate 
the LH 2 temperature and exit at 300K to avoid any ice on the set up. 

• A camera enables to see and record what happens into the cell, everything is 
recorded. 

• Flowmeters enable to control precisely the helium flow into the cryogenic circuit. 
However, they need to work at ambient temperature, that is why an heaters are 
placed before them , to set gas helium at room temperature. 

• A movable system that allows to check the LH 2 temperature along the vertical 
axis of the cell and to verify that one remains in isothermal conditions during the 
experiment is also installed. It is basically composed of rack and pinion. 

• A safety valve is also installed in the system. If the cell explodes, hydrogen is then 
released into the vacuum of the cryostat, leading to a strong pressure increase. 
Then the valve is automatically opened and the hydrogen is exhausted to the 
outside of the building. 

• A connection to the pumping system 
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Mechanical description of the cryostat 


The vacuum vessel of the cryostat is made of stainless steel 304L. It has been formed 
by boiler manufacturing. It is only 1 mm thick. That’s why there are several rein- 
forcements installed all over its entire height to prevent from crashing when making 
vacuum inside. Its volume is around 430 L and the medium diameter is 185 mm. 

A flat flange closes the cryostat. On the one hand, it mechanically supports all com- 
pound inside the cryostat as already mentioned and on the other hand it ensures 
airtightness between the inside and the outside of the cryostat. 

For electrical connection, Jaeger plugs are installed , enabling a simple an ergonomic 
access to the installation during maintenance. 

Both endoscopes can be manipulated independently one from the other. They are 
activated by electric motors controlled by a supervision panel on the computer. 

Mirrors are placed at the end of each endoscope. They are linked by cables which 
can be pulled by step by step motor to adjust their angular position. A spring is also 
installed in order to set the mirror back in position while keeping the cable free. 

Lighting endoscope has a removable frosted mirror, which enable to go from paral- 
lel to diffuse light. It can be activated by a switch on to of the cryostat. 

The diameter of endoscopes is D endo =A() mm and they are almost 2 m high. Hence, 
according to Eq (20) and (21) below, for a wavelength ranging from 400 and 800 nm, 
the smallest detail that can be seen is around 25-50 /rm 

3.1.2 The observation system 

The experimental cell is enclosed in the cryostat. Since everything is metallic, it is 
opaque, which makes the direct observation of the cell from the outside impossible. 
That is why an optical system is installed from the top flange. There is a light source 
(white light) whose the intensity can be modulated and a CCD camera with special lens 
to catch pictures. The observation is done by transmission of the light through the cell. 
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Endoscopes 


Since the camera and the lighting system are far away from the cell, an optical system 
had to be built for the light to follow a certain path, going from the lamp, through the 
cell and reaching the camera. An optic fiber guides the light path through the lighting 
endoscope until it reaches the cell. Before this, light rays going out of this fiber (and 
thus considered as punctual source) are going through several lenses for the light to 
diffuse in a parallel manner. 



Figure 16: Observation system 


CCD camera 


The camera is placed on the top flange of the vision endoscope to get the rays that 
are going out from the cell. A magnification objective is added for the rays to cover 
enough surface on the CCD sensor. 

The capture system is composed of : 

• an analogous and monochromatic CCD camera (Watec-902B) 

• a clock overlay (for date and time) 
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• a VCR (video cassette recorder) to VHS tape. 

• a printer top get an instant printout 

• a normal definition screen (SD) 


All compound are plugged in serial by RCA or BNC wires. The resolution is 768x494. 
The main advantage of this system is its quite simple installation. As soon as it is 
switched on, an image appears. However, post-processing of recordings are quite com- 
plex since it requires a numerical conversion. The poor resolution of recording doesn’t 
prevent from doing very sharp picture analysis. The camera is associated with a 70 
mm focal length objective with a focal doubler. Hence it’s a 140 mm focal length. 

Angular resolution 


Observation via an endoscope creates an alteration in quality of the real situation. 
The diffraction limits the angular resolution through the following formula: 

0«1,2— ^ — (20) 

Uendo 

where 9 is the maximal resolution power, \ (in m) the wavelength of the light, and 
D endo (in m), the diameter of the endoscope, corresponding to the numerical aperture. 
Since 9 is small, it leads to : 

tan9 » 9 sa — (21) 

Lendo 

where 1 (in m) is the size of the smallest element that can be seen by this system and 
L endo the length of the endoscope. 

3.1.3 The experimental cell 

There are several prerequisites to built an experimental cell for study of thermal transfer 
in microgravity: 

• It has to be transparent to be able to see what happens in there 

• It has to be regulated in temperature 

• It has to contain an heating component to apply an heating flux to the fluid 

• No leakage are allowed 
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That is the reason why the cell is built as presented on figure 17. A transparent sap- 
phire cylinder of 30 mm diameter is bridled at each side by one heat exchanger made 
of copper material. Heating resistances enable to regulate the temperature of the cell. 
Another advantage of the sapphire is its good conductivity at 20K. Its mechanical 
properties, such as tensile elastic limit, are compatible with the maximum possible 
pressure that are used. 

Stainless steel has very good mechanical properties, but poor thermal ones. It’s a mate- 
rial often used in cryogenics, unfortunately, its magnetic properties are disappointing. 
For our purpose, this very small interaction between material and magnetic field is 
totally advantageous. 


Heat exchangers are made of copper for better thermal conductivity. Other me- 
chanical pieces are made of stainless steel 304L, except for the inner ring presented 
below on figure 17 which is in titanium. 



Figure 17: Modelisation of the cell (on the left) and real cell (on the right) 

Titanium is chosen for this ring to thermally insulate the heating part from the 
bottom of the cell. It’s a ring 100 jum tick, shown in figure 18. The sealing of the cell 
is done with Indium seals. This has several advantages: on the one hand, it doesn’t 
become breakable at the temperature achieved by our experiment (20K) and it is inert 
with hydrogen. On the other hand, it very easy to shape it the needed way. as soon 
as it is installed, two clamps are squeezed. 
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Figure 18: section of the bottom of the cell (on the left) and picture of the bottom of 
the real cell (on the right) 

The cell itself (inside the sapphire) has a diameter of 30 mm and is 100 mm high. 
This cell was used for previous experiments with liquid oxygen. That’s why all copper 
elements have been recovered by a very thin layer (1/jrn) of gold. 

Regulation of the temperature 


Heat exchangers have been machined into the clamps (see fig 18). To carefully regulate 
the temperature, the technique utilized in our apparatus is to set a helium gas flow 
inside heat exchangers situated on both sides of the sapphire tube. As it has been 
experimentally observed, the temperature regulation is much better with gas helium 
than with liquid helium. 

The helium flow rate in each exchanger is imposed by the operator. As shown in figure 
13, the fluid arrives in a single pipe and is divided into the two ways near the cell. 
Heating resistance of about 100 Q are installed in each heat exchanger. Thermometers 
are constantly measuring the temperature of each of them. Their precision is about 
0.01K. The regulation is done by the Programmable Logic Controller (PLC) , which 
does a PID type regulation. 

Application of an heat flux 


At the center of the cell bottom flange, an heating resistor is installed inside the 
copper piece, enabling to heat it up only via Joule effect (see fig 18). Its exchange 
surface is 1 cm 2 . As explained above, this piece is thermally insulated from the rest 
of the exchanger with the titanium ring. The electrical generator gives 15 W to the 
resistor, which is sufficient since the theoretical heat flux is around 12 W.cm -2 [2] at 
normal temperature and pressure. 
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3.1.4 Thermal sensors 


To measure the temperature of the cell, there are several cryogenic CERNOX sensors. 
The heating element and each exchanger are equipped with such sensors. This ther- 
mometer is made of resistance whose value is proportional to the temperature. Hence, 
to determinate the temperature, it’s only necessary to measure the resistance and to 
associate the corresponding temperature value. Measures are 4-wire type: since errors 
due to resistance variation of the wire are reduced compared to a 2-wire type measure. 
The electrical connexion is then ensured by 4 wires with a diameter of 0.2 mm in man- 
ganin ( 86% copper, 12% manganese, and 2% nickel.). This material has been chosen 
thanks to its virtual zero temperature coefficient of resistance value and to minimize 
heat transfer between point that need electrical connections. 

Thermalisations are installed on the cold piece to limit even more the conduction ther- 
mal transfer. 

3.1.5 The magnetic map 

After having seen the material itself, it is interesting to detail the surrounding of the 
experiment, and especially the magnetic held resulting from the coil. 

A magnetic map has been established by using the software ’Levitation’ developed 
by the LNCMI. It uses a taylor expansion of the magnetic held to establish such a 
map. The hgure 19 below shows the residual gravity vector present in the considered 
volume, resulting from Eq (5) and Eq (6). 

The different curves correspond to the residual gravity, expressed in % of g. 
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The x-axis scale is 1 cm. Only half of the residual gravity map is presented above, 
since the second one is symmetric regarding to the coil axis. Half of the cell is repre- 
sented by the bold rectangle. Its bottom center is placed at the exact compensation 
point , at 12.6 cm above the center of the coil (0,0) 


3.2 Sizing of some elements 

In cryogenics, everything is sized. There is no place for chance. In this section the 
sizing of two different elements is presented: the current supply wires and the thermal 
heat exchanger. 
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3.2.1 Sizing of the current supply wire 


A 


Current supply wire 


L 



Cold piece that need to be 
regulated 


Resistor 



Figure 20: Drawing of the principle element needed to size current supply wire. 

Figure 20 presents the main element used to size current supply wire. A one dimension 
analysis is done. 

The power due to joule effect in a wire is 


with p the resistivity in il.rri, dx a small length element, I the current and A the section 
of the wire 

According to Eq(14), the power due to conduction in the wire is 


with A the thermal conductivity, A the section of the wire, and T the temperature, 
it leads to : 



(23) 


QdQ = p\I 2 dT 


thus: 
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with Q f the power brought at the cold part and Q c the power brought at the hot (the 
other) part of the wire. Moreover, Qf = Q CO nd + Q joules since in vacuum, all power due 
to joule effect is evacuated by conduction. The goal is to reduce heat input. Hence, 
Q c — 0 is imposed due to the thermalisation. Since we are looking for minimizing the 
heat flux at the cold part, we need Qcond = Q joule, hence: 

The Wiedemann-Frantz formula gives A p = LT, with A the thermal conductivity, p 
the resistivity T the temperature and L = 2.44 * 10 _8 VF.fhiF _2 a constant. It finally 
leads to : 

2 [ Tc A dT 

~A ~ yj L * (T c 2 - T7J ^ 

Since a thermalisation is done on a thermal shield at 80 K, the problem is splitted in 
two part: from ambient temperature at 300K to 80 K and from 80 to 20 K 


• For the first case: the needed power to bring is overestimated at 20 W. With a 
resistor of 100 fl, it corresponds to I =0.44 A. Tables give the conductivity integral 
between 300K and 80K : 76625,48 W.m~ l . Hence, we get : b = = 3.5 * 

10 -4 A/cm . The length of the wire L’ is 3m, then its diameter is : 



4> = 0.7 mm 


• The same reasoning gives a diameter of 0.2 mm between 80 and 20 K, where the 
length of the wire is 50 cm. 
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3.2.2 Sizing of the thermal shield heat exchanger 

A thermal shield is placed all around the cell in order to reduce the radiative heat 
coming from the hot (300K) vacuum vessel of the cryostat. 



Figure 21: Drawing of the principle elements in the vicinity of the thermal shield. 

As shown in figure 21, cold gas helium, coming from a tank ( not drawn on the right 
side of the figure), arrives into the cell to cool it down. Because of the heat exchanger 
pipe length, the helium temperature at the cell outlet is 20K. The helium gas is then 
utilized to cool down the thermal shield. In the following, one shows how this shield is 
size up. 

Moreover, an additional power of 2W is set on the cell during the experiment, to con- 
trol the temperature very precisely. 


• Calculation of the needed flow rate 


Firstly, the thermal shield exchanger is considered as perfect : the tempera- 
ture of the shield is equal to the temperature of the leaving helium. 

The heat balance of the cell has to be done : it receive radiating power from the 
thermal shield, the 2W additional power and the conduction power coming from 
the handling rod.: 


$1 — ^'2 a 6 screen ~ T Lll) + 2W + P cond 

with cr the Stephan constant, e the emissivity and S the area of the cell. For long 
cylinders, iq^ 2 is equal to: 
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r l ->-2 — 7 j 

62 + ~~ e ^ el 

For small emissivity, such as in our case (between 0.02 and 0.1 for copper and 
between 0.1 and 0.2 for steel), Fi_> 2 is often approximated by F\^ 2 ~ 0.5 
As this heat is evacuated with the heat exchanger, one gets : 


$ 2 = m'Cp(T 0 f - T if ) 


with m’ the helium flow rate, 2^ the temperature of the leaving helium, T t f the 
temperature of the incoming helium and C p the thermal capacity (in J/kg/K). 

Then, the heat balance of the thermal shield has to be done: it receives radiative 
power from the cryostat vacuum vessel. 

i 3 = F^aeS(Ti„ ng -Ti reen ) 

In our case, the previous relation becomes 

ts = 0.5aeS(Ti, lng -T^ en ) 


and evacuates it with its heat exchanger and by conduction into the handling rod: 


$ 4 = m'Cp(T ots - T of ) + P cond 


with T ots the temperature of the leaving helium. 

Pcond is calculated according to Eq(14). The values of the integral of conductivity, 
coming from tables has been linearised. 

Calculations have been done using an iterative method: first, the flowrate is 
chosen and imposed. Using (j ) i = (f) 2 the temperature of the thermal shield is 
calculated. 

This value is then put into the equation 0 3 = 0 4 . A new value of the flowrate is 
then obtained. This loop is done until reaching convergence. The final flowrate 
is 745 L/H of helium gas at normal temperature and pressure. 
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Calculation of the thermal shield exchanger temperature 


L 


T 4 Output He 



t 3 t 2 
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Ti Input He (20K) 


Figure 22: Drawing of the thermal shield exchanger: the core (in orange) and the pipe 
(in blue). 

The exchanged power is given by : 

05 — hSATi 0 g 

with h the convection coefficient, S the exchange area and 

\ rf-i (T2-T1)-(T3-T4) 

^J-log T2i T\ 


This power is also equal to 


05 = m'CpiTi - Ti) 


As same as previously, the method is iterative. The following system has to be 
solved: 


05 06 — 03 Pcond, 

By imposing a length of 359 mm, and the flowrate previously obtained, it corre- 
sponds to a final temperature of the thermal shield of 60K, using e = 0.1 for the 
copper emissivity 

3.3 Experimental protocol 

After having described the installation, its surrounding and the sizing of some element, 
let’s enter into the experimental part. A protocol is established to cool the installation 
down. It is described in the following chapter. 
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3.3.1 Starting the installation 

In this section, the protocol to prepare, run and stop an experiment of boiling crisis of 
liquid hydrogen in magnetic levitation is described. To understand this protocol, it is 
necessary to refer to different parts below on figure 23. 



The experimental protocol consists of the following steps: 

Simple checks: A number of checks are performed before all things: 

• Test all heaters, thermometers, flow meters and pumps of the cryostat. 
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Installation of the experimental cells: 


Before inserting the cell in the cryostat: 

• Check the good behaviour of thermometers and heating of the cell, 

• Check for any possible leaks of hydrogen and helium circuit of the cell by applying 
a vacuum in both circuits, then doing a leak detection with helium. The leakage 
rate must be less than 10 ~ 8 mb.l/s. Let the circuit at atmospheric pressure after 
the test. 

Installation of the cryostat inside the coil: 


The cryostat is inserted into the held hole of the coil with the overhead crane. This op- 
eration is very delicate. It is performed using the crane available in the hall of LNCMI. 
Since the cryostat is very high, movements of the crane can generate an important de- 
viation. The introduction of the cryostat must be very slow in order to keep everything 
safe inside. 

Vacuum in the cryostat: 


When the cryostat is placed in the hole , a secondary vacuum has to be proceeded, by 
connecting a pump to the valve PV100 (see figure 13). The vacuum in the cryostat 
must be of the order of 10 -5 mbar. The pressure gauge is given by the PT102. 

Rinsing and draining of the experimental cell: 


it is a necessity to purge the experimental cell under helium as follows: 

• Starting of the pump PP100 (pumping of the cell), 

• Opening the valves HV106 and HV103. Valves HV101, HV102, HV104 and 
HV105 are closed. The valve HV110 must always remain open, 

• When the pressure is low enough (PT101 = 10 _1 mb), the manual valve is operated 
closed HV106, 

• Opening of the helium bottle and valve HV101 until a pressure of 2 bar on PT101, 

• Closing of the valve HV101, 

• Repeat these steps 3 times. 

The cell must remain under helium. HV101 valve must remain open. The helium 
pressure is 2 bars. 
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Cooling and filling of the cell: 


• Connection of the helium siphon on the helium cooling circuit of the cell. 

• Regulation of the helium flow in the two exchangers of the cell to lower its tem- 
perature to approximately 20K. The average colling speed is 150 K/H. The cell 
is then filled as follows: 

- Closing of HV101 (HV103, HV106, HV107 and HV110 open, HV102, HV104 
and 105 closed HV), 

- Pumping the cell with the pump PP100, 

— Regulation of W101 on TT101 to get a temperature of 20K, 

- Regulation of W102 on TT102 to get a temperature of 20K, 

— Regulation of W103 on TT103 to get a temperature of 20K, 

- Closing of the valve PV100, 

- Opening of the bottle of hydrogen, 

— Opening of HV102 slowly so as to fill the cell with liquid hydrogen, 

- Adjustment to the desired pressure when the cell is full, . 

Endoscopes must be set, centered on the cell. 

3.3.2 Stopping the installation 

This last section is devoted to stopping the installation safely. The order of actions is 
crucial: 

• Closing of the valve HV102 

• Closing of the bottle of hydrogen. 

• Closing of HV100 and cut the flow of helium (FT201 and 202) 

• Pumping the hydrogen of the cell, 

• Stopping the regulations W101, W102, W103, W201 and W202 

• Putting the cell under helium after rinsing of the cell (see previous "Rinse and 
drain the experimental cell"), 

• Closing of the valve HV101 
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PS100 is a turbo-molecular pump, enabling to have a high secondary vacuum. It’s 
connected in serial with a primary pump. PS100 is very sensitive to magnetic field. 
Indeed, the rotor of the pump rotates at 36 000 RPM, which could create high eddy 
current leading to its destruction. That’s why it has been placed 6 meters away from 
the hole magnetic field. However, it creates much more pressure losses in the pumping 
pipes dur to their length. Hence the pumping of the cryostat is much slower. In 
this experiment it’s a dynamic vacuum: the pump is always running, even when the 
installation is hot (300K). Thus the pumping time is no longer a problem. 

3.3.3 Running of a measurement campaign 

After all these descriptions, it could be interesting here for the reader to remember the 
main objective of this campaign: the establishment of the Nukiyama curve for Liquid 
hydrogen ,(LH 2 ), (see figure 8) in microgravity. 

Finding of the exact compensation point: 


For experiments at normal gravity, as soon as the cell is at 20K, measurements can 
start. However, for microgravity of Og and O.lg, the exact compensation point has to 
be found. First of all, the cryostat is moved vertically by a lift table to place the cell at 
its theoretical position determined by the magnetic map (see figure 19). The current 
in the coil is risen until reaching its theoretical value. 

A small heat flux density is then applied to create some bublles. Figure 6 explained 
why hydrogen is radially unstable. Hence, at Og, bubbles should move radially. The 
observation of phenomenon is done directly by using the camera. 

If bubbles are still going up, the current is increased until getting the radial boil- 
ing. The speed of the bubble is also a precious information to find the compensation 
point: indeed, the closer to this point , the slower the bubbles. A scanning of a large 
band of current is done all around the point to be sure of having the real compensation 
point. 

In our case, this point was obtained for a current of 25100 A in the inner coil and 
20541 A in the outer one, corresponding to a magnetic field of 14.31 T at the magnetic 
center. The description of these coils is precisely done in Appendix section. 

The boiling curve is obtained by applying an heat flux density, it can be useful to 
remind that the overheating is defined by 

A T sat = T wall - T sat 

When the heat flux is imposed for the first time, the overheating can be very high 
before getting boiling: several nucleation sites can stay inactive. To activate all of 
then, a quite high flux , to the order of magnitude of the critical heat flux (CHF) is 
imposed. Measures can then start and are totally reproducible for a same heat flux. 
The temperature of the wall is measured by the thermal sensor TT100 (see it position 
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in figure 18 page 29), and the heat flux is imposed directly across the resistance, this 
power is entirely transmitted to the liquid. 

Whatever the gravity, measures and the detection of the CHF remain the same: it’s a 
thermal manner. 



Figure 24: Time evolution of the heating element (TT100) during the detection of CHF 

As shown on figure 24, the brutal increase of the temperature for an almost con- 
stant flux is characteristic of the CHF. 

As soon as the CHF is detected, the imposed flux is reduced to prevent the heating 
element to be destroyed. 

Procedure: 


This description will be done using figure 24. 

Let’s call T 0 a minimum temperature, smaller than the minimal temperature corre- 
sponding to the flux of film boiling mode. The heating compound is maintained to this 
temperature by applying a flux q 0 . 

At the time the flux is increased to q\ . the temperature evolves. At t < t 2 , 
when the temperature is stable (TT100=Ti), the point is taken. 

At the time t 2 , the flux is increased to q 2 . the temperature evolves. At t < t 2 , 
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when the temperature is stable (TT 100 =T 2 ), the point is taken. 


At the time t 3 , the flux is increased to q 3 . after a short time, the temperature in- 
creases rapidly: it’s the boiling crisis. T 3 is a value close to the maximum overheating. 
Hence, the critical flux is between : 


92 < qcHF < 93 

(25) 

the flux is then decreased to the first value. This procedure is ther 
a smaller value q 3 instead of q 3 . 

If the boiling crisis appears, then is obvious that 

i reproduced using 

92 < qcHF < q 3 

(26) 

Otherwise, 


9jt < Qchf < 9.3 

(27) 

The precision criterion has been chosen at 0.5W on the flux. 
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4 Results 


This master thesis has been done in collaboration with CEA Grenoble, LNCMI, Air 
liquide and CNES (the french Space Agency). Several measurements have been im- 
posed during experimental campaign, defined in a contract. There was a confidentiality 
clause in this contract, that’s why the curves are presented here without any scale on 
their axis, in order to keep confidentiality. 

Several Boiling curves had to be established for different thermodynamics properties 
that are presented on the following table: 


Level of gravity 

Pressure (bar) 

Temperature (K) 


1 

T sat = 20.23 


1 

T sat -2 = 18.23 

lg 

2 

T sat - 0.5 = 22.33 


2 

T sat -2 = 20.83 


1 

T sat -2= 18.23 

Og 

2 

T sat -2 = 20.83 


1 

T sat -2= 18.23 

0.1 g 

2 

T sat -2 = 20.83 


Table 1: Resume of the measurements that have to be performed 
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4.1 At lg 

4.1.1 lbar at saturation 

The first curve has been done for 1 bar and the corresponding saturation temperature. 



Figure 25: Nukiyama’s curve of liquid hydrogen at lbar and the corresponding satu- 
rated temperature under lg. 


The cell was half empty to be sure that the saturation condition were fulfilled. The 
pool boiling part of theoretical curve corresponds to Kutelatdze correlation [2]: 


r , kc^p}- 2 V - 75 A _ 

q = 1.8 * 10- 13 * f, , Ar 2 a : 

//0'625 (7 0.9 / \1.5pl.5 


(28) 


Where q is the heat flux (in W.cm~ 2 ), k is the thermal conductivity (in W.mT 2 .K~ l , 
C p is the specific heat capacity (in J.kg^.K* 1 ), p t is the density of the liquid, p is the 
pressure (in Pa), p is the viscosity in Pa.s, cr is the surface tension (in N/m), A is the 
latent heat of vaporisation (in J /kg) and p v is the density of the vapour. 

The film boiling part of the theoretical curve corresponds to Breen and Westwater 
correlation [2]: 


q = 


0.37 + 0.28( 

( 


.W* 

AT 


)1/S( ^f AT U/4 

gAp f } kjp f g\' 


(29) 


Where D is the characteristic length (in m). For this correlation, A' is the effective 
heat of vaporization given by [2]: 


, (A- 0.340C P AT) 2 

A “ A 


(30) 
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There have been two different sets of experiments which gave quite dispersive results. 

The cell has then be totally filled with liquid hydrogen. It enabled to impose a smaller 
temperature than the saturated one for a fixed pressure. Moreover, the temperature of 
the liquid pool is ever better controlled since both exchangers are in contact with the 
liquid. 

4.1.2 lbar and T=T sat -2K 

The following curve has been done with a temperature 2K smaller than the saturated 
one for lbar : T=18.23K 



Figure 26: Nukiyama’s curve of liquid hydrogen at lbar and T snr 2K under lg. 

As shown in figure 26 above, results are closer one from each other and also close 
to the theoretical curve. However it can be underlined that this theoretical one is the 
same as the previous one. Indeed, as it can be seen on Eq (28), since for liquid hydro- 
gen, the influence of the temperature is very small on different paramater (such as C p , 
etc...), the subcooling of this liquid is neglectible on the shape of the curve. 

There are some points for small overheating which are suspicious. They have been 
done at the very end of this set of experiment. Before emptying the cell, some small 
impurities have been seen at its bottom. It could have been some oxygen or nitrogen 
ice-cube present due to a small leak in the hydrogen feeding system. These impurities 
could be the reason of these last measurements. 

Then all the circuit has been controlled and all manual valves have been verified and 
tightened before starting the next set of experiments. 
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4.1.3 2 bars and T = T sat - 0.5K 


The pressure has then be increased to 2 bars and the temperature has been regulated 
to T=T sar 0.5 =22.33K. Results are shown below: 



Figure 27: Nukiyama’s curve of liquid hydrogen at 2bar and T sat -0.5K under lg. 

The theoretical curve has been established thanks to eq (28), with the data ex- 
tracted from NIST table. 

As it can be seen in figure 27 the experimental value of the CHF is much higher 
than the theoretical one. This measurements are also in general shifted to the left, i.e 
the overheating is smaller, the heat transfer is better. 
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4.1.4 2 bars and T = T sat - 2K 


The temperature has then be decreased to T = T sat - 2=20. 83K , keeping the same 
pressure. Results are shown below: 


.-V 


♦ Theoretical values 

♦ 1 set of experiment 

♦ 2nd set of experimen 


Figure 28: Nukiyama’s curve of liquid hydrogen at 2bar and T sat -2K under lg. 

As shown in figure 28, two sets of experiment have been carried out. Whereas both 
of them are similar in the him boiling part, in the nucleate boiling one, there are some 
discrepancies. These two sets have been done one day apart. 

In the second one the thermal transfer was much better since the overheating is shifted 
to the left. Moreover it can be seen that the CHF was found higher in the second 
one than in the first one, which is coherent with the previous affirmation. Indeed, the 
better the heat transfer, the higher the CHF. 
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4.2 At Og 

The magnetic field has then been set to obtain the microgravity conditions (see figure 
19 page 31). The exact compensation point of the gravity has been found experimen- 
tally by fixing the correct height and progressively increasing the current in the coil,i.e. 
increasing the magnetic field. The evolution of the shape of bubble are observed using 
the camera. 

As soon as bubbles grow radially and very slowly, the couple (I,z) is defined as the 
right one. Measurements can then start and results are shown in the following figures 

4.2.1 1 bar and T=T sat -2K 

The first set of measurements has been done for a pressure of 1 bar and a temperature 
2K smaller than its saturated one: T=18.23K . 



Figure 29: Nukiyama’s curve of liquid hydrogen at lbar and T sat - 2K under Og. 

The reproducibility of the results has been well obtained. Indeed as it can be seen 
in figure 29, for the same imposed heat flux, the overheating is almost exactly the same. 

The shape is as expected. However, even if in this master thesis no values can be 
expressed, it can be notified that the CHF is very different from the same configuration 
under normal gravity (lg). 
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Figure 30: Snapshots of the video of nucleate boiling under Og at lbar and T sat - 2K, 
from upper left corner to downer right one 

As it can be seen above, under Og, bubble are leaving radially the heater, to ag- 
glomerate on a bigger bubble at one side of the cell. 

figure 30 shows the him boiling in Og for these thermodynamic conditions: 



Figure 31: Snapshots of the video of nucleate boiling under Og at lbar and T sat - 2K 

As expected, the bubble is immobile. All the evaporated liquid is condensed at the 
interface. This big bubble is almost 10mm high and 20mm long. 
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4.2.2 2 bars and T=T sat -2K 


The pressure is then up to 2 bars and the temperature is regulated to T =Tsat — 
2=20. 83K. The following figure shows results. 



Figure 32: Nukiyama’s curve of liquid hydrogen at 2bar and T sat -2K under Og. 

In figure 32, it can be seen that here also the reproducibility is good. Which enable 
to assert that everythink is correct. Moreover, the expected shape is also obtained. 
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4.3 Under O.lg 

4.3.1 2 bars and T=T sat -2K 

To obtain O.lg, the current is calculated according to Eq (13), page 11. It led to an 
intensity of 23812 A in the inner coil and 19487 A in the outer coil. The resulting 
magnetic field at the center of both coils is 13.58 T. 



Figure 33: Nukiyama’s curve of liquid hydrogen at 2bar and T sat - 2K under O.lg. 

Obtained results are very satisfying, and the reproducibility is also been verified. 
The global shape fits totally as expected. 



Figure 34: Snapshots of the video of film boiling under O.lg at 2bar and T sar 2K 

In the first picture (upper left corner), the big leaving bubble can be seen and 
another one is forming. The size of the bubbles explicitly shows that we are in film 
boiling. Indeed, the bubble measure around 13 mm diameter ; which is higher than 
the heating compound size. The power injected by the heater is 4.9W 
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5 Discussion 

In the previous section, row have to be carried out to obtain real values for the heat 
flux and the measured temperature. 

5.1 Corrections 

5.1.1 Temperature of the heating part 



Figure 35: Position of the thermal sensor in the heating compound (green rectangle on 
the top). 

The real temperature of the heater is given by: 

AT rea i = T meas — A T cond 
AP e i 

with A T cond = ™ as c , and P meas = UI 

A c 7T L) heater 

All values are in SI units. Here e c i = 3.6 * 10 _3 m, A c — f(T meas ) and D heater = 
1.128* 10 -2 m. A T reat is the real temperature of the wall (in K). T rneas is the measured 
temperature of the sensor, AT cond is the temperature difference between the sensor and 
the wall, P m eas is the imposed power(in W) ,e cl is the distance between the sensor and 
the wall, Dheater is the diameter of the heating compound. 

This is the main correction that has to be applied to obtain true values tor the wall 
temperature of the heating area. 
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5.1.2 Pressure in the cell 

The real pressure in the cell is given by : 

Preal = Pmeas + AP, with A P = pn 2 -A.g.h 

Where pn 2 = 70.3 kg/m 3 at 1 bar and 20K, A is the percentage of g (0%, 10% or 100%), 
and h is the height of the cell (0.1m). 

The estimation of the height of the liquid hydrogen corresponds to the height of the 
sapphire tube, added with the height of the upper bridle. During experiment, this 
height can be higher. Hence, this correction Ap is the minimum one to apply to mea- 
sured pressures P meas - Corrections Ap are respectively 0.69, 0.069 and 0 mbar for lg, 
O.lg and Og. 

5.1.3 Power transfered to the liquid 

The real power transfered to the liquid via the heating area is (see figures 35 and 36 

Preal = Pmeas ~ W, - W 2 ~ W 3 (31) 

Where W\ (in W) is the losses in the titanium ring, W% (in W) is the losses in the disc 
and W% is the losses in the current supply wires. 



Figure 36: description of the different powers taking place in the heating compound. 


W\ = 


2Atjt7T G-ring {Preal Tbottom) 


ln{ 


Pring 
T ring 


(32) 
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Where A tit = /(^" — ^ bottom ) ^ 6ring = 0.1 * 10 3 m, R r i ng = 9-5 * 10 3 
frinq = 5.74 * 10“ 3 m. 


‘2\Gw'Ke diS c{T dit 


i / Rdisc \ 

'"W 


Where A G i 0 = f( Tdlsc ^ botto11x ), e disc — 5 * 10 3 m, R ring 
Tring = 1.55 * 10 -3 m. 


3.4 * 10" 2 


W 3 = - 


A c t xd 2 wire (T rea i + —ATcond — T therm ) 


= =2, l wire = 0.5 m, and d wire = 2 * 10 4 m, and 


A c — f{T reai + — A T cond — T therm ) 
e C 2 


It corresponds to the following maximal errors: 


Nucleate boiling 

Film boiling 

W lmax = 2.3 %P meas 
W 2max = 1.3 %P meas 
W 3ma x = 0.2 %P mea s 

Wi max = 10.4 %P meas 
W 2max = 1.05 %P meas 
W 3max = 0.4%P meas 


Table 2: Errors related to the measured power P mea s 


m and 


(33) 
m and 


(34) 
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5.2 Uncertainties on the losses related to the titanium ring 

Measurements of these uncertainties have been done while the cell was under vacuum. 
The temperature difference between each side of the ring (AT) is plotted as a function 
of the power. Figure 37 shows the difference between the theoretical curve and the 
obtained measurements. 



The maximum deviation from the theoretical curve correspond to 3% of the mea- 
sured value, which is totally acceptable. 
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Q (W/cm 2 ) 


5.3 Comparisons 

5.3.1 Effect of the gravity 

In the following condition of pressure and temperature, when the gravity level increases, 
it increases the critical heat flux and then the needed flux in film boiling for a given 
temperature is also higher : 



Figure 38: Comparison between different level of gravity for P =2 bar and T =T sat — 2 k 
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Q (W/cm 2 ) Q (W/cm 2 ) 


♦ Og 

_uJ9_ 


DT sat 

Figure 39: Comparison between different level of gravity for P =1 bar and T =T sat — 2k 

5.3.2 Effect of the pressure 

Under the following conditions of gravity and temperature, when the pressure increases, 
the critical heat flux also increases. In the film region, the effect of the pressure seems 
to be more significant for small gravity. 




Figure 40: Comparison between the different pressure with the same level of gravity 
(Og) and subcooling : T =T sat — 2 K 
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Q (W/cm 2 ) Q (W/cm 2 ) 



Figure 41: Comparison between the different pressure with the same level of gravity 
(lg) and subcooling : T =T sat — 2 K 

5.3.3 Effect of the subcooling 

In the following condition of pressure and gravity, when the subcooling increases, the 
critical heat flux doesn’t vary significantly. 



Figure 42: Comparison between the subcooling for the same pressure (2bar) and level 
of gravity (lg) 
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5.4 Problems encountered 

In this section, some problems that have been encountered during experiment are 
presented. 

5.4.1 Leakage in the cell feeding circuit 

There is a small capillary in which the hydrogen gas goes through before entering the 
cell. Before this entrance to the cell, there are roughly 12 m of pipe until the hydrogen 
bottle: 2.5m inside the cryostat and the rest of it outside. There are three handling 
valves in this circuit which, in case of defects, can lead to pollute the incoming hydrogen. 

It turned out that it happens for curve at lg. That is probably the reason why these 
results are very dispersive: The liquid hydrogen was not the same between the begin- 
ning and the end of the experiment, since small amount of air (oxygen+ nitrogen) was 
slowly diluted in it. 

Moreover, this leakage caused bigger problems: it plugged the capillary. Indeed, the 
solidification temperature of Oxygen is 60 K and the capillary was at 20K. It slowly 
created ice-plugs of oxygen. The speed of formation was certainly very slow due to 
the size of the leakage. Hence these ice-plugs were not present at the beginning of the 
experiment. 

The major problem appeared at the end of the set of experiment, when the cell would 
have been emptied : while increasing the temperature, the liquid became gas. How- 
ever, none of this gas could escape the cell, hence the pressure increased, following the 
phase diagram. The critical point for hydrogen is at 33K and the corresponding pres- 
sure is 13 bars. After this temperature, the fluid became supercritical. It was then 
impossible to exactly know what happened in the cell, even with the camera. 

A video of the disappearance of the interface between normal fluid and supercriti- 
cal fluid has been recorded, as it can be seen in figure 43. 

Finally the ice-plugs melt and was ejected at 50— 60K. Obviously, a huge peak of 
pressure has been measured while it happened. However, a security valve, regulated at 
3 bar, evacuated the hydrogen before damaging anything. 
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Figure 43: Disappearance of the interface at the critical point 

5.4.2 Fast breakdown of the magnetic field 

Because of the bad weather condition we had two brutal breakdowns of the current 
(and thus the magnetic held, going from 14T to OT in 2s). 

It created eddy currents in the copper exchangers of the cell and in the copper thermal 
shield all around. It resulted in a huge amount of energy (700 000 J), which instanta- 
neously evaporated the whole cell (totally filled with liquid hydrogen). The pressure 
increased at such a point that the indium seal crept, letting the hydrogen gas going 
into the vacuum of the cryostat. Lorentz’ force that resulted in this breakdown inflated 
the thermal shield and made it in its plastic stress. Then to remove this shield, an 
hydraulic jack has been used, as shown in the picture 44. 200 bars were applied to the 
piston, corresponding to a force of 14 300 N. 
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Figure 44: Use of an hydraulic jack to extract the thermal shield from the cryostat. 

5.4.3 Losses in the cooling circuit 

The previous problem damaged the shield. While a new one was made, the previ- 
ous one was cut at the bottom and adapted to be able to continue experiments. An 
aluminium part was put to reduce radiation coming from the cryostat vacuum vessel 
(300K) as much as possible. However, It totally changed the behaviour of the cooling 
system, which was sized for the previous one. 

The helium flow rate became insufficient to cool the cell. The understanding of what 
happened was no straightforward. The first idea has been to evaluate the losses in the 
gas siphon. This has been done by imposing a flow rate with a flowmeter placed di- 
rectly at the output of the siphon. The temperature of the flow was measured. Hence, 
by using the formula: 

$ = m'CpiTout - Tin) 

Where $ is the power, m’ the flowrate and T in and T out the temperature at the en- 
trance (in the helium tank, ie: 4.2K) and at the out of the siphon. Theses losses were 
estimated to 7 W, which gave a second explanation of the limitation in the cooling. 

To be able to pursue our experiments, a pumping system has been installed at the 
output of the cooling (helium) system. Hence, the applied (relative) pressure in the 
pipes went from 0.5 bar to 1.5 bar, which enabled to increase the flowrate in the circuit 
and to continue cooling without any limitation. 
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6 Conclusion 


This master thesis enabled to have a good understanding of boiling heat transfers in 
microgravity. Different thermodynamic parameters have been studied independently 
one from each other to be able to draw solid conclusion regarding the behaviour of 
the fluid (liquid hydrogen) in such conditions. It also enabled to apply lots of different 
physics : mechanical behaviour at 20 K, cryogenics, vacuum technique, computational 
physics, magnetostatic behaviour, heat transfers... 

It turns out that at Og, the CHF is much lower than at 0.1 or lg, whereas this CHF 
follows the increase of pressure. 

Results are very satisfying and the main goal is achieved: This study can be used 
to design the new type of motor for Ariane spacecraft, that can be turned on into 
space after a while of inactivity. 
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